High-resolution HLA genotyping improves PIRCHE-II assessment of molecular mismatching in kidney transplantation.

HLA matching in solid organ transplant is performed with the aim of assessing immunologic compatibility in order to avoid hyperacute rejection and assess the risk of future rejection events. Molecular mismatch algorithms are intended to improve granularity in histocompatibility assessment and risk stratification. PIRCHE-II uses HLA genotyping to predict indirectly presented mismatched donor HLA peptides, though most clinical validation studies rely on imputing high resolution (HR) genotypes from low resolution (LR) typing data. We hypothesized that use of bona fide HR typing could overcome limitations in imputation, improving accuracy and predictive ability for donor-specific antibody development and acute rejection. We performed a retrospective analysis of adult and pediatric kidney transplant donor/recipient pairs (N = 419) with HR typing and compared the use of imputed LR genotyping verses HR genotyping for PIRCHE-II analysis and outcomes. Imputation success was highly dependent on the reference population used, as using historic Caucasian reference populations resulted in 10 % of pairs with unsuccessful imputation while multiethnic reference populations improved successful imputation with only 1 % unable to be imputed. Comparing PIRCHE-II analysis with HR and LR genotyping produced notably different results, with 20 % of patients discrepantly classified to immunologic risk groups. These data emphasize the importance of using multiethnic reference panels when performing imputation and indicate HR HLA genotyping has clinically meaningful benefit for PIRCHE-II analysis compared to imputed LR typing.

Integrative multi-omics profiling in human decedents receiving pig heart xenografts.

In a previous study, heart xenografts from 10-gene-edited pigs transplanted into two human decedents did not show evidence of acute-onset cellular- or antibody-mediated rejection. Here, to better understand the detailed molecular landscape following xenotransplantation, we carried out bulk and single-cell transcriptomics, lipidomics, proteomics and metabolomics on blood samples obtained from the transplanted decedents every 6 h, as well as histological and transcriptomic tissue profiling. We observed substantial early immune responses in peripheral blood mononuclear cells and xenograft tissue obtained from decedent 1 (male), associated with downstream T cell and natural killer cell activity. Longitudinal analyses indicated the presence of ischemia reperfusion injury, exacerbated by inadequate immunosuppression of T cells, consistent with previous findings of perioperative cardiac xenograft dysfunction in pig-to-nonhuman primate studies. Moreover, at 42 h after transplantation, substantial alterations in cellular metabolism and liver-damage pathways occurred, correlating with profound organ-wide physiological dysfunction. By contrast, relatively minor changes in RNA, protein, lipid and metabolism profiles were observed in decedent 2 (female) as compared to decedent 1. Overall, these multi-omics analyses delineate distinct responses to cardiac xenotransplantation in the two human decedents and reveal new insights into early molecular and immune responses after xenotransplantation. These findings may aid in the development of targeted therapeutic approaches to limit ischemia reperfusion injury-related phenotypes and improve outcomes.

Exploration of potential immune mechanisms in cervical dystonia.

BACKGROUND: Although there are many possible causes for cervical dystonia (CD), a specific etiology cannot be identified in most cases. Prior studies have suggested a relationship between autoimmune disease and some cases of CD, pointing to possible immunological mechanisms. OBJECTIVE: The goal was to explore the potential role of multiple different immunological mechanisms in CD. METHODS: First, a broad screening test compared neuronal antibodies in controls and CD. Second, unbiased blood plasma proteomics provided a broad screen for potential biologic differences between controls and CD. Third, a multiplex immunoassay compared 37 markers associated with immunological processes in controls and CD. Fourth, relative immune cell frequencies were investigated in blood samples of controls and CD. Finally, sequencing studies investigated the association of HLA DQB1 and DRB1 alleles in controls versus CD. RESULTS: Screens for anti-neuronal antibodies did not reveal any obvious abnormalities. Plasma proteomics pointed towards certain abnormalities of immune mechanisms, and the multiplex assay pointed more specifically towards abnormalities in T lymphocytes. Abnormal immune cell frequencies were identified for some CD cases, and these cases clustered together as a potential subgroup. Studies of HLA alleles indicated a possible association between CD and DRB1*15:03, which is reported to mediate the penetrance of autoimmune disorders. CONCLUSIONS: Altogether, the association of CD with multiple different blood-based immune measures point to abnormalities in cell-mediated immunity that may play a pathogenic role for a subgroup of individuals with CD.

Belatacept with time-limited tacrolimus coimmunosuppression modifies the 3-year risk of eplet mismatch in kidney transplantation.

Solid organ transplant donor-recipient eplet mismatch has been correlated with donor-specific antibody (DSA) formation, antibody-mediated rejection, and overall rejection rates. However, studies have been predominantly in patients on tacrolimus-based immunosuppression regimens and have not fully explored differences in ethnically and racially diverse populations. Evidence indicates that patients on belatacept have lower rates of DSA formation, suggesting mediation of the immunogenicity of mismatched human leukocyte antigen polymorphisms. We performed a retrospective, single-center analysis of class II eplet disparity in a cohort of kidney transplant recipients treated using belatacept with tacrolimus induction (Bela/TacTL) or tacrolimus regimens between 2016 and 2019. Bela/TacTL (n = 294) and tacrolimus (n = 294) cohorts were propensity score-matched with standardized difference <0.15. Single-molecule eplet risk level was associated with immune event rates for both groups. In Cox regression analysis stratified by eplet risk level, Bela/TacTL immunosuppression was associated with a decreased rate of DSA (hazard ratio [HR] = 0.4), antibody-mediated rejection (HR = 0.2), and rejection (HR = 0.45). In the low-risk group, cumulative graft failure was lower for patients on Bela/TacTL (P < .02). Analysis of eplet mismatch burden may be a useful adjunct in identifying high-risk populations with increased immunosuppression requirements and should encourage the design of allocation rules to incentivize lower-risk pairings without negatively impacting equity in access.

HLA sensitization is associated with an increased risk of primary graft dysfunction after heart transplantation.

Primary graft dysfunction (PGD) is a leading cause of early morbidity and mortality following heart transplantation (HT). We sought to determine the association between pretransplant human leukocyte antigen (HLA) sensitization, as measured using the calculated panel reactive antibody (cPRA) value, and the risk of PGD. METHODS: Consecutive adult HT recipients (n = 596) from 1/2015 to 12/2019 at 2 US centers were included. Severity of PGD was based on the 2014 International Society for Heart and Lung Transplantation consensus statement. For each recipient, unacceptable HLA antigens were obtained and locus-specific cPRA (cPRA-LS) and pre-HT donor-specific antibodies (DSA) were assessed. RESULTS: Univariable logistic modeling showed that peak cPRA-LS for all loci and HLA-A was associated with increased severity of PGD as an ordinal variable (all loci: OR 1.78, 95% CI: 1.01-1.14, p = 0.025, HLA-A: OR 1.14, 95% CI: 1.03-1.26, p = 0.011). Multivariable analysis showed peak cPRA-LS for HLA-A, recipient beta-blocker use, total ischemic time, donor age, prior cardiac surgery, and United Network for Organ Sharing status 1 or 2 were associated with increased severity of PGD. The presence of DSA to HLA-B was associated with trend toward increased risk of mild-to-moderate PGD (OR 2.56, 95% CI: 0.99-6.63, p = 0.053), but DSA to other HLA loci was not associated with PGD. CONCLUSIONS: Sensitization for all HLA loci, and specifically HLA-A, is associated with an increased severity of PGD. These factors should be included in pre-HT risk stratification to minimize the risk of PGD.

IgM marks persistent IgG anti-human leukocyte antigen antibodies in highly sensitized heart transplant patients.

BACKGROUND: Sensitization to human leukocyte antigens (HLA) is a persistent problem in heart transplant (HT) candidates. We sought to characterize the anti-HLA antibody and circulating B cell repertoire in a cohort of highly sensitized HT candidates. METHODS: We assessed immunoglobulin G (IgG) and immunoglobulin M (IgM) anti-HLA antibodies using Luminex single antigen bead assays in a cohort of 11 highly sensitized (HS; calculated panel reactive antibody ≥ 90%) and 3 mildly sensitized (MS) candidates. We also performed B cell receptor repertoire sequencing (BCRseq) in HS candidates and 33 non-candidate controls. HLA antibody strength was measured by mean fluorescence intensity (MFI). RESULTS: We found that IgM anti-HLA antibodies were present in all HS candidates, but with a lower breadth and strength as compared to IgG. When anti-HLA IgG specificities intersected with IgM, binding strength was higher. In contrast, there were IgM but no intersecting IgG specificities for the MS group. In four candidates in the HS group, IgG anti-HLA antibodies decreased in both breadth and strength after HT, but the decrease in strength was smaller if the IgG possessed a specificity that intersected with pre-transplant IgM. BCRseq revealed larger B cell clonotypes in HS candidates but similar diversity as compared to controls. CONCLUSIONS: IgM marks IgG anti-HLA antibodies with higher strength before HT and persistence after HT. The presence of IgM intersecting IgG for an anti-HLA specificity may be a useful approach to determine which donor HLA should be avoided for a sensitized candidate.

Combined imputation of HLA genotype and self-identified race leads to better donor-recipient matching.

Allogeneic Hematopoietic Cell Transplantation (HCT) is a curative therapy for hematologic disorders and often requires human leukocyte antigen (HLA)-matched donors. Donor registries have recruited donors utilizing evolving technologies of HLA genotyping methods. This necessitates in-silico ambiguity resolution and statistical imputation based on haplotype frequencies estimated from donor data stratified by self-identified race and ethnicity (SIRE). However, SIRE has limited genetic validity and presents a challenge for individuals with unknown or mixed SIRE. We present MR-GRIMM "Multi-Race Graph IMputation and Matching" that simultaneously imputes the race/ethnic category and HLA genotype using a SIRE based prior. Additionally, we propose a novel method to impute HLA typing inconsistent with current haplotype frequencies. The performance of MR-GRIMM was validated using a dataset of 170,000 donor-recipient pairs. MR-GRIMM has an average 20 % lower matching error (1-AUC) than single-race imputation. The recall metric (sensitivity) of the race/ethnic category imputation from HLA was measured by comparing the imputed donor race with the donor-provided SIRE. Accuracies of 0.74 and 0.55 were obtained for the prediction of 5 broad and 21 detailed US population groups respectively. The operational implementation of this algorithm in a registry search could help improve match predictions and access to HLA-matched donors.

Genetic variation in the human leukocyte antigen region confers susceptibility to Clostridioides difficile infection.

Clostridioides difficile (C. diff.) infection (CDI) is a leading cause of hospital acquired diarrhea in North America and Europe and a major cause of morbidity and mortality. Known risk factors do not fully explain CDI susceptibility, and genetic susceptibility is suggested by the fact that some patients with colons that are colonized with C. diff. do not develop any infection while others develop severe or recurrent infections. To identify common genetic variants associated with CDI, we performed a genome-wide association analysis in 19,861 participants (1349 cases; 18,512 controls) from the Electronic Medical Records and Genomics (eMERGE) Network. Using logistic regression, we found strong evidence for genetic variation in the DRB locus of the MHC (HLA) II region that predisposes individuals to CDI (P > 1.0 × 10-14; OR 1.56). Altered transcriptional regulation in the HLA region may play a role in conferring susceptibility to this opportunistic enteric pathogen.

Genotype List String 1.1: Extending the Genotype List String grammar for describing HLA and Killer-cell Immunoglobulin-like Receptor genotypes.

The Genotype List (GL) String grammar for reporting HLA and Killer-cell Immunoglobulin-like Receptor (KIR) genotypes in a text string was described in 2013. Since this initial description, GL Strings have been used to describe HLA and KIR genotypes for more than 40 million subjects, allowing these data to be recorded, stored and transmitted in an easily parsed, text-based format. After a decade of working with HLA and KIR data in GL String format, with advances in HLA and KIR genotyping technologies that have fostered the generation of full-gene sequence data, the need for an extension of the GL String system has become clear. Here, we introduce the new GL String delimiter "?," which addresses the need to describe ambiguity in assigning a gene sequence to gene paralogs. GL Strings that do not include a "?" delimiter continue to be interpreted as originally described. This extension represents version 1.1 of the GL String grammar.

Molecular HLA mismatching for prediction of primary humoral alloimmunity and graft function deterioration in paediatric kidney transplantation.

Introduction: Rejection remains the main cause of allograft failure in paediatric kidney transplantation and is driven by donor-recipient HLA mismatching. Modern computational algorithms enable assessment of HLA mismatch immunogenicity at the molecular level (molecular-mismatch, molMM). Whilst molMM has been shown to correlate with alloimmune outcomes, evidence demonstrating improved prediction performance against traditional antigen mismatching (antMM) is lacking. Methods: We analysed 177 patients from the CERTAIN registry (median follow-up 4.5 years). molMM scores included Amino-Acid-Mismatch-Score (AAMS), Electrostatic-Mismatch-Score (EMS3D) and netMHCIIpan (netMHC1k: peptide binding affinity ≤1000 nM; netMHC: binding affinity ≤500 nM plus rank <2%). We stratified patients into high/low-risk groups based on risk models of DSA development. Results: Donor-specific HLA antibodies (DSA) predominantly targeted the highest scoring molMM donor antigen within each HLA locus. MolMM scores offered superior discrimination versus antMM in predicting de novo DSA for all HLA loci; the EMS3D algorithm had particularly consistent performance (area under the receiver operating characteristic curve (AUC) >0.7 for all HLA loci vs. 0.52-0.70 for antMM). ABMR (but not TCMR) was associated with HLA-DQ molMM scores (AAMS, EMS3D and netMHC). Patients with high-risk HLA-DQ molMM had increased risk of graft function deterioration (50% reduction in baseline eGFR (eGFR50), adjusted HR: 3.5, 95% CI 1.6-8.2 high vs. low EMS3D). Multivariable modelling of the eGFR50 outcome using EMS3D HLA-DQ stratification showed better discrimination (AUC EMS3D vs. antMM at 2 years: 0.81 vs. 0.77, at 4.5 years: 0.72 vs. 0.64) and stratified more patients into the low-risk group, compared to traditional antMM. Conclusion: Molecular mismatching was superior to antigen mismatching in predicting humoral alloimmunity. Molecular HLA-DQ mismatching appears to be a significant prognostic factor for graft function deterioration in paediatric kidney transplantation.

GRAMM: A new method for analysis of HLA in families.

Recently, haplo-identical transplantation with multiple HLA mismatches has become a viable option for stem cell transplants. Haplotype sharing detection requires the imputation of donor and recipient. We show that even in high-resolution typing when all alleles are known, there is a 15% error rate in haplotype phasing, and even more in low-resolution typings. Similarly, in related donors, the parents' haplotypes should be imputed to determine what haplotype each child inherited. We propose graph-based family imputation (GRAMM) to phase alleles in family pedigree HLA typing data, and in mother-cord blood unit pairs. We show that GRAMM has practically no phasing errors when pedigree data are available. We apply GRAMM to simulations with different typing resolutions as well as paired cord-mother typings, and show very high phasing accuracy, and improved allele imputation accuracy. We use GRAMM to detect recombination events and show that the rate of falsely detected recombination events (false-positive rate) in simulations is very low. We then apply recombination detection to typed families to estimate the recombination rate in Israeli and Australian population datasets. The estimated recombination rate has an upper bound of 10%-20% per family (1%-4% per individual).

Corrigendum: Molecular HLA mismatching for prediction of primary humoral alloimmunity and graft function deterioration in paediatric kidney transplantation.

[This corrects the article DOI: 10.3389/fimmu.2023.1092335.].

Machine learning reveals limited contribution of trans-only encoded variants to the HLA-DQ immunopeptidome.

Human leukocyte antigen (HLA) class II antigen presentation is key for controlling and triggering T cell immune responses. HLA-DQ molecules, which are believed to play a major role in autoimmune diseases, are heterodimers that can be formed as both cis and trans variants depending on whether the α- and ß-chains are encoded on the same (cis) or opposite (trans) chromosomes. So far, limited progress has been made for predicting HLA-DQ antigen presentation. In addition, the contribution of trans-only variants (i.e. variants not observed in the population as cis) in shaping the HLA-DQ immunopeptidome remains largely unresolved. Here, we seek to address these issues by integrating state-of-the-art immunoinformatics data mining models with large volumes of high-quality HLA-DQ specific mass spectrometry immunopeptidomics data. The analysis demonstrates highly improved predictive power and molecular coverage for models trained including these novel HLA-DQ data. More importantly, investigating the role of trans-only HLA-DQ variants reveals a limited to no contribution to the overall HLA-DQ immunopeptidome. In conclusion, this study furthers our understanding of HLA-DQ specificities and casts light on the relative role of cis versus trans-only HLA-DQ variants in the HLA class II antigen presentation space. The developed method, NetMHCIIpan-4.2, is available at https://services.healthtech.dtu.dk/services/NetMHCIIpan-4.2 .

HLA amino acid Mismatch-Based risk stratification of kidney allograft failure using a novel Machine learning algorithm.

OBJECTIVE: While associations between HLA antigen-level mismatches (Ag-MM) and kidney allograft failure are well established, HLA amino acid-level mismatches (AA-MM) have been less explored. Ag-MM fails to consider the substantial variability in the number of MMs at polymorphic amino acid (AA) sites within any given Ag-MM category, which may conceal variable impact on allorecognition. In this study we aim to develop a novel Feature Inclusion Bin Evolver for Risk Stratification (FIBERS) and apply it to automatically discover bins of HLA amino acid mismatches that stratify donor-recipient pairs into low versus high graft survival risk groups. METHODS: Using data from the Scientific Registry of Transplant Recipients, we applied FIBERS on a multiethnic population of 166,574 kidney transplants between 2000 and 2017. FIBERS was applied (1) across all HLA-A, B, C, DRB1, and DQB1 locus AA-MMs with comparison to 0-ABDR Ag-MM risk stratification, (2) on AA-MMs within each HLA locus individually, and (3) using cross validation to evaluate FIBERS generalizability. The predictive power of graft failure risk stratification was evaluated while adjusting for donor/recipient characteristics and HLA-A, B, C, DRB1, and DQB1 Ag-MMs as covariates. RESULTS: FIBERS's best-performing bin (on AA-MMs across all loci) added significant predictive power (hazard ratio = 1.10, Bonferroni adj. p < 0.001) in stratifying graft failure risk (where low-risk is defined as zero AA-MMs and high-risk is one or more AA-MMs) even after adjusting for Ag-MMs and donor/recipient covariates. The best bin also categorized more than twice as many patients to the low-risk category, compared to traditional 0-ABDR Ag mismatching (∼24.4% vs âˆ¼ 9.1%). When HLA loci were binned individually, the bin for DRB1 exhibited the strongest risk stratification; relative to zero AA-MM, one or more MMs in the bin yielded HR = 1.11, p < 0.005 in a fully adjusted Cox model. AA-MMs at HLA-DRB1 peptide contact sites contributed most to incremental risk of graft failure. Additionally, FIBERS points to possible risk associated with HLA-DQB1 AA-MMs at positions that determine specificity of peptide anchor residues and HLA-DQ heterodimer stability. CONCLUSION: FIBERS's performance suggests potential for discovery of HLA immunogenetics-based risk stratification of kidney graft failure that outperforms traditional assessment.

HLA Class I Genotype Is Associated with Relapse Risk after Allogeneic Stem Cell Transplantation for NPM1-Mutated Acute Myeloid Leukemia.

Mutation-bearing peptide ligands from mutated nucleophosmin-1 (NPM1) protein have been empirically found to be presented by HLA class I in acute myeloid leukemia (AML). We hypothesized that HLA genotype may impact allogeneic hematopoietic stem cell transplantation (allo-HCT) outcomes in NPM1-mutated AML owing to differences in antigen presentation. We evaluated the effect of the variable of predicted strong binding to mutated NPM1 peptides using HLA class I genotypes from matched donor-recipient pairs on transplant recipients' overall survival (OS) and disease-free survival (DFS) as part of the primary objectives and cumulative incidence of relapse and nonrelapse mortality (NRM) as part of secondary objectives. Baseline and outcome data reported to the Center for International Blood and Marrow Transplant Research from a study cohort of adult patients (n = 1020) with NPM1-mutated de novo AML in first (71%) or second (29%) complete remission undergoing 8/8 matched related (18%) or matched unrelated (82%) allo-HCT were analyzed retrospectively. Class I alleles from donor-recipient pairs were analyzed for predicted strong HLA binding to mutated NPM1 using netMHCpan 4.0. A total of 429 (42%) donor-recipient pairs were classified as having predicted strong-binding HLA alleles (SBHAs) to mutated NPM1. In multivariable analyses adjusting for clinical covariates, the presence of predicted SBHAs was associated with a lower risk of relapse (hazard ratio [HR], .72; 95% confidence interval [CI], .55 to .94; P = .015). OS (HR, .81; 95% CI, .67 to .98; P = .028) and DFS (HR, .84; 95% CI, .69 to 1.01; P = .070) showed a suggestion of better outcomes if predicted SBHAs were present but did not meet the prespecified P value of <.025. NRM did not differ (HR, 1.04; P = .740). These hypothesis-generating data support further exploration of HLA genotype-neoantigen interactions in the allo-HCT context.

Unrelated Stem Cell Donor HLA Match Likelihood in the US Registry Incorporating HLA-DPB1 Permissive Mismatching.

Donor-recipient HLA matching at the DPB1 locus improves the outcomes of hematopoietic stem cell transplantation (HCT). Retrospective outcome studies found that in HCTs matched for all 8 alleles of the A, B, C, and DRB1 loci at high resolution (8/8 match), few transplantations were also allele-matched at the DPB1 locus. DPB1 allele matching was once thought to be logistically impractical; however, a DPB1-permissive mismatch model based on T cell epitope (TCE) reactivity expands the proportion of suitable donors. To understand the likelihood of finding a DPB1-permissive donor, we sought to expand population genetic match likelihood models for the US unrelated donor registry, the National Marrow Donor Program (NMDP). After extending HLA haplotype frequency estimates to include the DPB1 locus, our models found that the likelihood of having a DPB1-permissive donor was not much lower than likelihood of 8/8 matching. A maximum of 5 additional donors would need to be typed to find a more optimal DPB1-permissive donor at least 90% of the time. Linkage disequilibrium patterns between the DPB1 locus and other classical HLA loci varied markedly by haplotype and population, indicating that the known recombination hotspot between DQ and DP gene complexes has not had a uniform impact; thus, DPB1-permissive donors are easier to identify within minority populations. DPB1 TCE categories were highly predictable from HLA typing at other loci when imputed with extended haplotype frequency data. Our overall results indicate that registry search strategies that seek a more optimally matched HCT donor encompassing HLA-DPB1 permissibility are likely to be highly productive.

Incidence of Pediatric Celiac Disease Varies by Region.

INTRODUCTION: The Environmental Determinants of Diabetes in the Young study follows an HLA risk selected birth cohort for celiac disease (CD) development using a uniform protocol. Children under investigation come from 6 different regions within Europe and the United States. Our aim was to identify regional differences in CD autoimmunity and CD cumulative incidence for children born between 2004 and 2010. METHODS: Children (n = 6,628) with DQ2.5 and/or DQ8.1 were enrolled prospectively from birth in Georgia, Washington, Colorado, Finland, Germany, and Sweden. Children underwent periodic study screening for tissue transglutaminase antibodies and then CD evaluation per clinical care. Population-specific estimates were calculated by weighting the study-specific cumulative incidence with the population-specific haplogenotype frequencies obtained from large stem cell registries from each site. RESULTS: Individual haplogenotype risks for CD autoimmunity and CD varied by region and affected the cumulative incidence within that region. The CD incidence by age 10 years was highest in Swedish children at 3%. Within the United States, the incidence by age 10 years in Colorado was 2.4%. In the model adjusted for HLA, sex, and family history, Colorado children had a 2.5-fold higher risk of CD compared to Washington. Likewise, Swedish children had a 1.4-fold and 1.8-fold higher risk of CD compared with those in Finland and Germany, respectively. DISCUSSION: There is high regional variability in cumulative incidence of CD, which suggests differential environmental, genetic, and epigenetic influences even within the United States. The overall high incidence warrants a low threshold for screening and further research on region-specific CD triggers.

Divergent HLA variations and heterogeneous expression but recurrent HLA loss-of- heterozygosity and common HLA-B and TAP transcriptional silencing across advanced pediatric solid cancers.

The human leukocyte antigen (HLA) system is a major factor controlling cancer immunosurveillance and response to immunotherapy, yet its status in pediatric cancers remains fragmentary. We determined high-confidence HLA genotypes in 576 children, adolescents and young adults with recurrent/refractory solid tumors from the MOSCATO-01 and MAPPYACTS trials, using normal and tumor whole exome and RNA sequencing data and benchmarked algorithms. There was no evidence for narrowed HLA allelic diversity but discordant homozygosity and allele frequencies across tumor types and subtypes, such as in embryonal and alveolar rhabdomyosarcoma, neuroblastoma MYCN and 11q subtypes, and high-grade glioma, and several alleles may represent protective or susceptibility factors to specific pediatric solid cancers. There was a paucity of somatic mutations in HLA and antigen processing and presentation (APP) genes in most tumors, except in cases with mismatch repair deficiency or genetic instability. The prevalence of loss-of-heterozygosity (LOH) ranged from 5.9 to 7.7% in HLA class I and 8.0 to 16.7% in HLA class II genes, but was widely increased in osteosarcoma and glioblastoma (~15-25%), and for DRB1-DQA1-DQB1 in Ewing sarcoma (~23-28%) and low-grade glioma (~33-50%). HLA class I and HLA-DR antigen expression was assessed in 194 tumors and 44 patient-derived xenografts (PDXs) by immunochemistry, and class I and APP transcript levels quantified in PDXs by RT-qPCR. We confirmed that HLA class I antigen expression is heterogeneous in advanced pediatric solid tumors, with class I loss commonly associated with the transcriptional downregulation of HLA-B and transporter associated with antigen processing (TAP) genes, whereas class II antigen expression is scarce on tumor cells and occurs on immune infiltrating cells. Patients with tumors expressing sufficient HLA class I and TAP levels such as some glioma, osteosarcoma, Ewing sarcoma and non-rhabdomyosarcoma soft-tissue sarcoma cases may more likely benefit from T cell-based approaches, whereas strategies to upregulate HLA expression, to expand the immunopeptidome, and to target TAP-independent epitopes or possibly LOH might provide novel therapeutic opportunities in others. The consequences of HLA class II expression by immune cells remain to be established. Immunogenetic profiling should be implemented in routine to inform immunotherapy trials for precision medicine of pediatric cancers.

Pathogenicity and impact of HLA class I alleles in aplastic anemia patients of different ethnicities.

Acquired aplastic anemia (AA) is caused by autoreactive T cell-mediated destruction of early hematopoietic cells. Somatic loss of human leukocyte antigen (HLA) class I alleles was identified as a mechanism of immune escape in surviving hematopoietic cells of some patients with AA. However, pathogenicity, structural characteristics, and clinical impact of specific HLA alleles in AA remain poorly understood. Here, we evaluated somatic HLA loss in 505 patients with AA from 2 multi-institutional cohorts. Using a combination of HLA mutation frequencies, peptide-binding structures, and association with AA in an independent cohort of 6,323 patients from the National Marrow Donor Program, we identified 19 AA risk alleles and 12 non-risk alleles and established a potentially novel AA HLA pathogenicity stratification. Our results define pathogenicity for the majority of common HLA-A/B alleles across diverse populations. Our study demonstrates that HLA alleles confer different risks of developing AA, but once AA develops, specific alleles are not associated with response to immunosuppression or transplant outcomes. However, higher pathogenicity alleles, particularly HLA-B*14:02, are associated with higher rates of clonal evolution in adult patients with AA. Our study provides insights into the immune pathogenesis of AA, opening the door to future autoantigen identification and improved understanding of clonal evolution in AA.